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Many energetic 7-ray pulsars discovered by Fermi are promising candidates for optical follow-ups. 
We present the results of the first deep optical observations of the two Vela-like Fermi pulsars PSR 
J1357-6429 and PSR J1048-5832 performed with the VLT. However, they have not been detected 
down to V~27 and V~27.6, respectively (3ct). These upper limits suggest an efficiency in converting 
spin-down power into optical luminosity < 7 x 10 -7 and < 6 x 10 -6 , respectively, lower than the 
Crab pulsar and, possibly, more compatible with the spin-down age of these two pulsars. 



INTRODUCTION 



Optical observations of pulsars are crucial to study 
the intrinsic properties of neutron stars, from the 
structure and composition of the interior, to the prop- 
erties and geometry of the magnetosphere. Histor- 
ically, X-ray and 7-ray observations have paved the 
way to the pulsar optical identifications, with all the 
presently identified pulsars also detected at high en- 
ergies (Mignani 2011). In particular, 5 out of the 7 7- 
ray pulsars detected by NASAs CGRO satellite have 
also been detected in the optical, suggesting that 7- 
ray detections highlight promising candidates for op- 
tical observations, since the emission at both ener- 
gies seems to correlate with the strength of the mag- 
netic field at the light cylinder (e.g., Shearer & Golden 
2001). The launch of the Fermi Gamma-ray Space 
Telescope opened new perspectives with the detec- 
tion of over 80 7-ray pulsars (Abdo et al. 2011). For 
most Fermi pulsars, however, no deep optical obser- 
vations have been carried out so far. In the North- 
ern hemisphere, an exploratory survey with 2.5m/4m 
telescopes at the La Palma Observatory has been per- 
formed (Shearer et al., these proceedings), while in 
the Southern hemisphere a survey with the VLT is in 
progress (see Mignani et al. 2011a). Here, we present 
the search for the optical counterparts of two Fermi 
pulsars (Mignani et al. 2011b). PSR J1357-6429 is 
a very young (7.31 kyr old) pulsar (P = 166.1 ms), 
at a distance of ~2.5 kpc, with a rotational energy 



loss rate E ~ 3 x 10 Jb erg s 



It has been de- 



tected as a 7-ray pulsar both by AGILE (Pellizzoni et 
al. 2009) and Fermi (Lemoine-Goumard et al. 2011), 
while its associated pulsar wind nebula (PWN) has 
been detected at TeV energies by HESS (Abramowski 
et al. 2011). PSR J1048-5832 is a 20.3 kyr old, 
Vela-like, pulsars (P=123.6 ms), tentatively detected 
by the CGRO (Kaspi et al. 2000) and now confirmed 
by Fermi (Abdo et al. 2009), at a distance of ^2.7 
kpc and with a£~ 2 x 10 36 erg s _1 comparable to 
that of PSR J1357-6429. Both pulsars are also de- 
tected in X-rays by both XMM-Newton and Chandra 
(Esposito et al. 2007; Lemoine-Goumard et al. 2011; 
Chang et al. 2011; Gonzalez et al. 2006) . No deep op- 
tical observations of these two pulsars have been ever 
reported. We used multi-band VLT images available 
in the ESO public archive to search for their optical 
emission 



II. OBSERVATIONS 

Optical images of the PSR J1357-6429 and PSR 
J1048-5832 fields were obtained with the VLT be- 
tween April 2009 and February 2010. Observations 
were performed in service mode with FORS2, a multi- 
mode camera for imaging and long-slit/multi-object 
spectroscopy (MOS). FORS2 was equipped with its 
red-sensitive MIT detector, a mosaic of two 2k?4k 
CCDs optimised for wavelengths longer than 6000 A. 
In its standard resolution mode, the detector has a 
pixel size of 0.25 arcsec (2x2 binning) which cor- 
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FIG. 1: VLT/F0RS2 30x30 arcsec image cutouts of PSR J1357-6429 (left) and PSR J1048-5832 (right) taken through 
the V filter. North to the top, East to the left, in size. The circles marks the pulsar radio (T=timing; I=interferometry) 
and Chandra positions. 



responds to a fieldofview of 8.3 x 8.3 arcmin over 
the CCD mosaic. Observations were performed with 
the standard low gain, fast read-out mode and in 
high-resolution mode (0.125 arcsec/pixel) for PSR 
J1357— 6429 (V,R,I bands) and in standard resolu- 
tion mode (0.25" /pixel) for PSR J1048-5832 (V 
band). Bright stars close to the PSR J1357-6429 
position were masked using the MOS slitlets. To al- 
low for cosmic ray removal and minimise saturation 
of bright stars in the field, sequences of short expo- 
sures (from 200 to 750 s) were obtained per each tar- 
get and per each filter. The total integration time 
was 14700 s (V), 8800 s (R), and 1600 s (I) for PSR 
J1357-6429 and 24000 s (V) for PSR J1048-5832. 
Exposures were taken in dark time and photomet- 
ric conditions, close to the zenith (airmass<1.3), 
and sub-arcsec image quality. We reduced the data 
through standard packages in IRAF for bias subtrac- 
tion, and flat-field correction using the closest-in-time 
bias and twilight flat-fields frames available in the 
ESO archive. Per each band, we aligned and average- 
stacked the reduced science images using the IRAF 
task drizzle applying a 3a filter on the single pixel 
average to filter out residual hot and cold pixels and 
cosmic ray hits. We applied the photometric cali- 
bration by using the extinction-corrected night zero 
points computed by the FORS2 pipeline and available 
through the instrument data quality control database 
(www.eso.org/qc). To register the pulsar positions on 
the FORS2 frames as precisely as possible, we re- 
computed their astrometric solution. We measured 
the star centroids through Gaussian fitting using the 
GAIA tool and used the code ASTR0M to compute the 
pixcl-to-sky coordinate transformation. For both pul- 
sar fields, we estimate that the overall (la) uncer- 



tainty of our FORS2 astrometry is ~ 0.2 arcsec. 



III. RESULTS 

As a reference to compute the pulsar positions on 
the FORS2 frames we used their most recently pub- 
lished radio and Chandra coordinates (see Mignani 
et al. 2011b). The radio and Chandra positions of 
PSR J1357-6429 and PSR J1048-5832 are shown in 
Fig. 1. For the former the Chandra position, averaged 
over multi-epoch consistent measurements, differs by 
1.2 ± 0.4 arcsec from the radio intcrfcrometry one, an 
amount which can be explained by the ^7.2 yr time 
span between the epochs for a pulsar transverse veloc- 
ity of 2100 ±700 km s _1 . For the latter, the difference 
between the Chandra and the radio timing positions 
amounts to 4.1 ±1.3 arcsec and cannot be explained by 
the pulsar proper motion between the two epochs (<~ 
5.4 yrs apart). However, the Chandra position is con- 
sistent with the radio-interferometry one, suggesting 
that the radio timing position was affected by timing 
noise. Since the multi-epoch Chandra coordinates ap- 
pear more reliable, in the following we assume them as 
a reference, although we conservatively evaluate pos- 
sible candidates at the radio positions. 

The radio position of PSR J1357-6429 (Fig. 1, 
left) falls close (0.65 arcsec) to a relatively bright field 
object (object A; V=21.8±0.06). We verified whether 
object A can be considered a possible counterpart to 
the pulsar. The chance coincidence probability which 
turns out to be ~0.01, i.e. not statistically compelling 
yet. If it were its counterpart, the flux of object A 
would imply an optical emission efficiency for PSR 
J1357— 6429 up to 10 times larger than the Crab pul- 
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sar and up to ^5600 times larger than the Vela pul- 
sar, to which it is more similar in its spin parame- 
ters. Thus, it is unlikely that object A is the pul- 
sar counterpart. We also compared object As colours 
(V=21.80±0.06; R=20.65±0.03; I=19.75±0.03) with 
those of field stars and we found that it has no peculiar 
colours which may suggest its association with the pul- 
sar. For instance, for a flat power-law spectrum, like 
that of the Vela pulsar, we would expect a (V-R)~0.2 
and a (VT)~0.1. Thus, we conclude that object A is 
unrelated to PSR J1357— 6429. The same arguments 
also rule out object B (V=23.05±0.13; R=21.86±0.07; 
I=20.5±0.05) as a candidate counterpart. We found 
a marginal evidence of a flux enhancement over the 
background at the edge of the Chandra error circle, 
which might come from the presence of a faint source 
(I « 24.6). However, the local crowding, with two 
relatively bright stars within a radius of ~ 1.5 arc- 
sec from the Chandra position, and the lack of detec- 
tions in the V and R-bands, make it problematic to 
determine whether such an enhancement come from a 
background fluctuation, perhaps produced by the PSF 
wings of the two stars, or it is associated with a real 
source and hence to a putative pulsar counterpart. 



The radio-timing position of PSR J1048-5832 (Fig. 
1, right) falls close to a V^24 object (C). Again, the 
chance coincidence probability is ^0.04, i.e. certainly 
not low enough to statistically claim an association. 
This is also ruled out by the pulsars energetics. The 
flux of object C would imply an optical emission effi- 
ciency for PSR J1048— 5832 up to 60 times larger than 
the Crab pulsar and up to ~200000 times larger than 
the Vela pulsar. Thus, it is extremely unlikely that 
object C is associated with PSR J1048-5832. Then, 
we searched for a possible counterpart at the pulsar 
Chandra position. No object is detected within, or 
close to, the Chandra (0.55 arcsec) error circle apart 
from a faint object (Star D; V ~ 26.7) visible ~ 1.6 
arcsec southwest of it. However, the offset is about 
three times the la uncertainty on the pulsar position. 
Thus, we deem the association unlikely both on the 
basis of the loose positional coincidence and on statis- 
tical grounds, with a chance coincidence probability 
of P ~ 0.08. Moreover, the lack of colour information 
makes it impossible to constrain the object's nature. 



No other possible counterparts are detected for both 
PSR J1357-6429 and PSR J1048-5832. We con- 
clude that they are unidentified in the VLT images 
and we computed the 3a upper limits on their flux. 
For PSR J1357-6429 we derive V>27.3, R>26.8, 
and I>24.6, while for PSR J1048-5832 we derived 
V>28. Moreover, we could not find evidence of coun- 
terparts of the X-ray PWNe detected by Chandra (see, 
Mignani et al. 2011b). 



IV. DISCUSSION 



We compared our V-band 1 flux upper limits with 
the pulsars' rotational energy loss rates. For PSR 
J1357— 6429, our upper limit of V <~ 27 corresponds 
to an optical luminosity upper limit L opt ~ 0.6- 
21.6 x 10 29 erg s" 1 , for a distance d = 2.4 ± 0.6 kpc 
and for an interstellar extinction Ay = 2.2~^\''[, based 
upon the Nh derived from the fit to the X-ray spec- 
tra.. This implies an emission efficiency upper limit 
Vopt ~ 0.2-7 x 10~ 7 . This value is at least a factor 
of 5 lower than the Crab pulsar and, possibly, closer 
to that of the Vela pulsar. For PSR J1048-5832 our 
upper limit of V ~ 27.6 implies (for d = 2.7±0.35 kpc 
and Ay = 5± 2 ; 2 ) upper limits ofL opt ~ 0.4-12.5 x 10 30 
erg s _1 and rj opt ~ 1.8-62.5 x 10~ 7 . In principle, this 
does not rule out optical emission efficiency compa- 
rable to that of the Crab pulsar, although the pulsar 
spin-down age (20.3 kyr) might also suggest in this 
case, Vela-like emission efficiency. 

We also compared the flux upper limits with the ex- 
trapolations in the optical domain of the X and 7-ray 
spectra of PSR J1357-6429 and PSR J1048-5832 . 
For PSR J1357— 6429, we assumed an X-ray power 
law (PL) with photon index Tx = 1-4 ± 0.5 plus a 
blackbody (BB) with temperature kT = 0.16+^ keV 
{N H = 0.4±°; 3 x 10 22 cm" 2 ; Esposito et al. 2007), and 
a 7-ray PL with photon index T 7 = 1.54 ± 0.41, and 
exponential cut-off at ~ 0.8 GeV (Lemoine-Goumard 
et al. (2011), while for PSR J1048-5832, we as- 
sumed a PL with photon index Fx = 2.4 ± 0.5 
(N H = 0.9t° ,t x 10 22 cm" 2 ; Marelli et al. (2011), and 
a 7-ray PL with photon index T 7 = 1.38 ± 0.13 and 
exponential cut-off at - 2.3 GeV (Abdo et al. 2009). 

The multi-wavelength spectral energy distributions 
(SEDs) of the two pulsars are shown in Fig. 2. For 
PSR J1357-6429 (Fig. 2, left) we see that the optical 
flux upper limits can be compatible with the extrap- 
olation of the X-ray power-law (PL). The optical flux 
upper limits are also compatible, with the possible ex- 
ception of the R-band one, with the extrapolation of 
the 7-ray PL, which does not allow us to prove that 
there is a break in the optical/7-ray spectrum, as ob- 
served in other pulsars. The multi-wavelength SED 
is different in the case of PSR J1048-5832 (Fig. 2, 
right), for which the optical V-band upper limit is 
compatible with the extrapolation of the steep X-ray 
PL only for the highest values of the Njj , but is well 
above the extrapolation of the flat 7-ray one. Interest- 
ingly enough, at variance with PSR J1357— 6429, no 
single model can describe the optical-to-7-ray mag- 
netospheric emission of PSR J1048-5832. Thus, the 
comparison between the multi-wavelength SEDs of 
these two Vela-like pulsars suggests that the occur- 
rence of spectral breaks might not correlate with the 
pulsar's age. 
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FIG. 2: Upper limits on the extinction-corrected optical fluxes of PSR J1357-6429 (left) and PSR J1048-5832 (right) 
compared with the low-energy extrapolations of the X-ray (solid black) and 7-ray (solid blue) spectral models that 
best fit the XMM-Newton, Chandra, and Fermi data. The dotted lines in the left panel correspond to the PL and BB 
components to the model X-ray spectrum. In both panels, the yellow and blue-shaded areas indicate the la uncertainty 
on the extrapolations of the X and 7-ray PLs, respectively. The vertical red dashed lines mark the uncertainty on the 
extinction-corrected optical flux upper limits, computed around the best-fit value of the Nh- 



V. SUMMARY AND CONCLUSIONS 

We used archival VLT/ FORS2 observations to 
perform the first deep optical investigations of 
the two Fermi pulsars PSR J1357-6429 and PSR 
J1048— 5832. We re-assessed the positions of the two 
pulsars from the analyses of all the available Chan- 
dra observations and a comparison with the published 
radio positions. For PSR J1357— 6429, this yielded 
a tentative proper motion of \i = 0.17 ± 0.055 arc- 
sec yr _1 (70° ± 15° position angle), which needs to 
be confirmed by future radio-interferometry observa- 
tions. For both pulsars, none of the objects detected 
around the Chandra positions can be considered vi- 
able candidate counterparts on the basis of their rel- 
atively large optical flux, > 3er offset from the pulsar 
position, and lack of peculiar colours with respect to 
the field stars. We then determined 3tr upper lim- 
its of V ~ 27 and V ~ 27.6 for PSR J1357-6429 
and PSR J1048— 5832, respectively, corresponding 
to optical emission efficiencies f] pt < 7 x 10" 7 and 
Vopt ~ 6xl0~ 6 . The VLT observations presented here 



are close to the sensitivity limits achievable with 10m- 
class telescopes under sub-arcsec seeing conditions. In 
the case of PSR J1357— 6429, deep, high spatial res- 
olution images with the HST are probably required 
to firmly claim an object detection from the flux en- 
hancement seen at the Chandra position (see section 
3). In the case of PSR J1048-5832, the large inter- 
stellar extinction (Ay ~ 5) hamper observations in 
the optical/near- ultraviolet. Deep observations in the 
near-infrared, either with the HST or with adaptive 
optic device at lOm-class telescopes, might represent 
a better opportunity to spot a candidate counterpart 
to the pulsar. 
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